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Abstract 

A comprehensive review will be given about the rich mathemati- 
cal structure of mean field spin glass theory, mostly developed, until 
now, in the frame of the methods of theoretical physics, based on deep 
physical intuition and hints coming from numerical simulation. Cen- 
tral to our treatment is a very simple and yet powerful interpolation 
method, allowing to compare different probabilistic schemes, by us- 
ing convexity and positivity arguments. In this way we can prove the 
existence of the thermodynamic limit for the free energy density of 
the system, a long standing open problem. Moreover, in the frame 
of a generalized variational principle, we can show the emergency of 
the Derrida-Ruelle random probability cascades, leading to the form 
of free energy given by the celebrated Parisi Ansatz. All these re- 
sults seem to be in full agreement with the mechanism of spontaneous 
replica symmetry breaking as developed by Giorgio Parisi. 



* e-mail: francesco.guerra@romal.infn.it 
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1 Introduction 



The mean field model for spin glasses, introduced by David Sherrington and 
Scott Kirkpatrick more that thirty years ago [I], 0, is a celebrated model. 
Hundreds and hundreds of articles have been devoted to its study during the 
years, appearing in the theoretical physics literature. 

The relevance of the model stems surely from the fact that it is intended 
to represent some important features of the physical spin glass systems, of 
great interest for their peculiar properties, exhibiting a new magnetic phase, 
where magnetic moments are frozen into disordered equilibrium orientations, 
without any long-range order. See for example [3] for a very readable review 
about the physical properties of spin glasses. 

But another important source of interest is connected with the fact that 
disordered systems, of the Sherrington-Kirkpatrick type, and their general- 
izations, seems to play a very important role for theoretical and practical 
assessments about hard optimization problems, as it is shown for example 
by Mark Mezard, Giorgio Parisi and Riccardo Zecchina in 0j. 

It is interesting to remark that the original paper was entitled "Solvable 
Model of a Spin-Glass", while a previous draft, as I have it from David 
Sherrington, contained even the stronger denomination "Exactly Solvable". 
However, it turned out that the very natural solution devised by the authors 
was valid only at high temperatures, or large external magnetic fields. While, 
at low temperatures, the proposed solution exhibited a nonphysical drawback 
given by a negative entropy, as properly recognized by the authors in their 
very first paper. 

It took some years to find an acceptable solution. This was done by 
Giorgio Parisi in a series of papers, by marking a radical departure from 
the previous methods. In fact, a very deep method of "spontaneous replica 
symmetry breaking" was developed. As a consequence the physical content 
of the theory was encoded in a functional order parameter of new type, and 
a remarkable structure emerged for the pure states of the theory, a kind of 
hierarchical, ultrametric organization. These very interesting developments, 
due to Giorgio Parisi, and his coworkers, are explained in a lucid way in the 
classical book Part of this structure will be recalled in the following. 

It is important to remark that Parisi solution is presented in the form 
of an ingenious and clever Ansatz. Until few years ago it was not known 
whether this Ansatz would give the true solution for the model, in the so 
called thermodynamic limit, when the size of the system becomes infinite, or 
it would be only a very good approximation for the true solution. 

The general structure offered by the Parisi solution, and their possible 
generalizations for similar models, exhibit an extremely rich and interesting 
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mathematical content. Very appropriately, Michel Talagrand has inserted a 
strongly suggestive sentence in the title to his recent book : "Spin glasses: 
a challenge for mathematicians" . 

As a matter of fact, how to face this challenge is a very difficult problem. 
Here we would like to recall the main features of a very powerful method, yet 
extremely simple in its very essence, based on a comparison and interpolation 
argument on sets of Gaussian random variables. 

The method found its first simple application in 0, where it was shown 
that the Sherrington-Kirkpatrick replica symmetric approximate solution was 
a rigourous lower bound for the quenched free energy of the system, uni- 
formly in the size. Then, it was possible to reach a long waited result jH]: the 
convergence of the free energy density in the thermodynamic limit, by an in- 
termediate step where the quenched free energy was shown to be subadditive 
in the size of the system. 

Moreover, still by interpolation on families of Gaussian random variables, 
the first mentioned result was extended to give a rigorous proof that the 
expression given by the Parisi Ansatz is also a lower bound for the quenched 
free energy of the system, uniformly in the size [§]. The method gives not 
only the bound, but also the explicit form of the correction in a quite in- 
volved form. As a recent and very important result, along the task of facing 
the challenge, Michel Talagrand has been able to dominate these correction 
terms, showing that they vanish in the thermodynamic limit. This milestone 
achievement was firstly announced in a short note [TO] , containing only a 
synthetic sketch of the proof, and then presented with all details in a long 
paper to be published on Annals of Mathematics ^1] . 

The interpolation method is also at the basis of the far reaching gener- 
alized variational principle proven by Michel Aizenman, Robert Sims and 
Shannon Starr in [T2j . 

In our presentation, we will try to be as self-contained as possible. We will 
give all definitions, explain the basic structure of the interpolation method, 
and show how some of the results are obtained. We will concentrate mostly 
on questions connected with the free energy, its properties of subadditivity, 
the existence of the infinite volume limit, and the replica bounds. 

For the sake of comparison, and in order to provide a kind of warm up, we 
will recall also some features of the standard elementary mean field model of 
ferromagnetism, the so called Curie- Weiss model. We will concentrate also 
here on the free energy, and systematically exploit elementary comparison 
and interpolation arguments. This will show the strict analogy between the 
treatment of the ferromagnetic model and the developments in the mean field 
spin glass case. Basic roles will be played in the two cases, but with different 
expressions, by positivity and convexity properties. 
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The organization of the paper is as follows. In Section 2, we introduce the 
ferromagnetic model and discuss behavior and properties of the free energy 
in the thermodynamic limit, by emphasing, in this very elementary case, 
the comparison and interpolation methods that will be also exploited, in a 
different context, in the spin glass case. 

Section 3 is devoted to the basic features of the mean field spin glass 
models, by introducing all necessary definitions. 

In Section 4, we introduce the Gaussian comparison and interpolation 
method, by giving simple applications to the existence of the infinite volume 
limit of the quenched free energy and to the proof of general variational 
bounds, by following the useful strategy developed in [T2"] . 

Section 5 will briefly recall the main features of the Parisi representation, 
and will state the main theorem concerning the free energy. 

Finally, Section 6 will be devoted to conclusions and outlook for future 
developments. 

It is a pleasure to thank the Organizing Committee, and in particular 
Professor Ari Laptev, for the kind invitation to talk in a such stimulating 
cultural athmosphere. 

2 The mean field ferromagnetic model. Struc- 
ture and results. 

The mean field ferromagnetic model is among the simplest models of statis- 
tical mechanics. However, it contains very interesting features, in particular 
a phase transition, characterized by spontaneous magnetization, at low tem- 
peratures. We refer to standard textbooks, for example [T3], for a full treat- 
ment, and a complete appreciation of the model in the frame of the theory 
of ferromagnetism. Here we consider only some properties of the free energy, 
easily obtained through comparison methods. 

The generic configuration of the mean field ferromagnetic model is defined 
through Ising spin variables Cj = ±1, attached to each site i = 1, 2, . . . , N. 

The Hamiltonian of the model, in some external field of strength h, is 
given by the mean field expression 



(m) * 

Here, the first sum extends to all N(N — l)/2 site couples, an the second to 
all sites. 
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For a given inverse temperature /?, let us now introduce the partition 
function Z N (f3, h) and the free energy per site Jn(P, h), according to the well 
known definitions 

Z N (J3,h)= exp{-/3H N (<r,h)), (2) 

(T1...(TJV 

-(3f N (f3,h) =N~ 1 E\ogZ N (f3,h). (3) 
It is also convenient to define the average spin magnetization 

™=^E^- ( 4 ) 

i 

Then, it is immediately seen that the Hamiltonian in JIJ can be equiva- 
lently written as 

H N (a,h) = -^Nm 2 -hJ2^ (5) 

i 

where an unessential constant term has been neglected. In fact we have 

5> i0 i = \ E a ^ = \ N2m2 ~ \ N > ( 6 ) 

{id) ij'&i 

where the sum over all couples has been equivalently written as one half the 
sum over all i,j with % ^ j, and the diagonal terms with i = j have been 
added and subtracted out. Notice that they give a constant because of = 1. 

Therefore, the partition function in (J2J) can be equivalently substituted 
by the expression 

Z N (f3,h)= exp(-^iVm 2 )exp(/3/ i ^a i ), (7) 

CTl. ..<JN i 

which will be our starting point. 

Our interest will be in the liniAr^oo iV -1 log Z N (fi, h). To this purpose, let 
us establish the important subadditivity property, holding for the splitting of 
the big N site system in two smaller N\ site and A^ 2 site systems, respectively, 
with N = Ni + N 2 , 

log Z N (P, h) < log Z Nl (P, h) + log Z N2 (0, h). (8) 

The proof is very simple. Let us denote, in the most natural way, by 
C7i, ... , cttvj the spin variables for the first subsystem, and by 0^+1, • • • , <?n 
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the N 2 spin variables of the second subsystem. Introduce also the subsys- 
tem magnetizations mi and m 2 , by adapting the definition (j3J) to the smaller 
systems, in such a way that 

Nm = N x mi + N 2 m 2 . (9) 

Therefore, we see that the large system magnetization m is the linear convex 
combination of the smaller system ones, according to the obvious 

m = -jt™<i + jy m 2- ( 10 ) 

Since the mapping m — > m 2 is convex, we have also the general bound, 
holding for all values of the a variables 

m < —m 1 + j^m 2 . (11) 

Then, it is enough to substitute the inequality in the definition (JZj) of Zn(/3, h), 
and recognize that we achieve factorization with respect to the two subsys- 
tems, and therefore the inequality Z N < Z Nl Z N2 . So we have established ©• 
From subadditivity, the existence of the limit follows by a simple argument, 
as explained for example in In fact, we have 



lim N~ l log Z N ((3, h) =mfN~ 1 fog Z N ((3,h). (12) 

N^oo N 

Now we will calculate explicitely this limit, by introducing an order parameter 
M, a trial function, and an appropriate variational scheme. In order to 
get a lower bound, we start from the elementary inequality m 2 > 2mM — 
M 2 , holding for any value of m and M. By inserting the inequality in the 
definition (|7J) we arrive at a factorization of the sum over cr's. The sum 
can be explicitely calculated, and we arrive immediately to the lower bound, 
uniform in the size of the system, 

N' 1 log Z N (/3,h) > log2 + logcosh/3(/i + M) - i/3M 2 , (13) 

holding for any value of the trial order parameter M. Clearly it is convenient 
to take the supremum over M. Then we establish the optimal uniform lower 
bound 

N- 1 log Z N (l3,h) > sup(log2 + log cosh /?(/i + M) - \l3M 2 ). (14) 
M 2 

It is simple to realize that the supremum coincides with the limit as 
N — > oo. To this purpose we follow the following simple procedure. Let 
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us consider all possible values of the variable m. There are N + 1 of them, 
corresponding to any number K of possible spin flips, starting from a given 
cr configuration, K = 0,1, . . . , N. Let us consider the trivial decomposition 
of the identity, holding for any m, 

1 = y^mM, (15) 

M 

where M in the sum runs over the N + 1 possible values of m, and 5 is 
Kroneker delta, beeing equal to 1 if M = N, and zero otherwise. Let us now 
insert (fTK|) in the definition (JZj) of the partition function inside the sum over 
cr's, and invert the two sums. Because of the forcing m = M given by the 
5, we can write m 2 = 2mM — M 2 inside the sum. Then if we neglect the 5, 
by using the trivial 5 < 1, we have un upper bound, where the sum over cr's 
can be explicitily performed as before. Then it is enough to take the upper 
bound with respect to M, and consider that there are N + 1 terms in the 
now trivial sum over M, in order to arrive at the upper bound 

N~ x log Z N {(3,h) < sup(log2 + logcosh/5(/i + M)-^/3M 2 ) + A^- 1 log(iV + l). 
M 2 

(16) 

Therefore, by going to the limit as iV — > oo, we can collect all our results 
in the form of the following theorem giving the full characterization of the 
thermodynamic limit of the free energy. 

Theorem 1. For the mean field ferromagnetic model we have 

lim N~ l log Z N (J3, h) = inf N' 1 log Z N (J3, h) (17) 

N^oo N 

= sup(log2 + logcosh/3(/i + M) — -f3M 2 ). (18) 
M 2 

This ends our discussion about the free energy in the ferromagnetic model. 
Now we are ready to attack the much more difficult spin glass model. But it 
will be surprising to see that, by following a simple extension of the methods 
here described, we will arrive to similar results. 

3 The basic definitions for the mean field spin 
glass model 

As in the ferromagnetic case, the generic configuration of the mean field spin 
glass model is defined through Ising spin variables Oi = ±1, attached to each 
site i = 1,2, . . . ,N. 
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But now there is an external quenched disorder given by the N(N — l)/2 
independent and identical distributed random variables J^, defined for each 
couple of sites. For the sake of simplicity, we assume each Jjy to be a centered 
unit Gaussian with averages E(Jij) = 0, E(Jfj) = 1. By quenched disorder 
we mean that the J have a kind of stochastic external influence on the system, 
without partecipating to the thermal equilibrium. 

Now the Hamiltonian of the model, in some external field of strength h, 
is given by the mean field expression 

H N (a, h, J) = — -== ^2 ■ , ij (T ' r! i - h^cTi. (19) 

Here, the first sum extends to all site couples, an the second to all sites. 
Notice the necessary to ensure a good thermodynamic behavior to the 
free energy. 

For a given inverse temperature (3, let us now introduce the disorder 
dependent partition function Z N (f3,h, J) and the quenched average of the 
free energy per site /jv(/3, h), according to the definitions 

Z N ((3,h,J)= exp(-(3H N (a,h,J)), (20) 

a 1 ...a N 

-PMP, h) = N~ l E log Z N ((3, h, J). (21) 

Notice that in (pTTf the average E with respect to the external noise is made 
after the log is taken. This procedure is called quenched averaging. It 
represents the physical idea that the external noise does not partecipate to 
the thermal equilibrium. Only the ex's are thermalized. 

For the sake of simplicity, it is also convenient to write the partition 
function in the following equivalent form. First of all let us introduce a family 
of centered Gaussian random variables k(o), indexed by the configurations 
a, and characterized by the covariances 

E(K(a)K(a')) = q 2 (a,a'), (22) 

where q(a, a 1 ) are the overlaps between two generic configurations, defined 
by 

q(a,cr') = N- 1 J2^ (23) 

i 

with the obvious bounds —1 < q(cr, a') < 1, and the normalization q(a, a) = 
1. Then, starting from the definition (|19|l . it is immediately seen that the par- 
tition function in (|2T)|) can be also written, by neglecting unessential constant 
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terms, in the form 

Z N (f3,h,J)= exp(P^ K (a))eMPhJ2°i), (24) 
which will be the starting point of our treatment. 

4 Gaussian comparison and applications 

Our basic comparison argument will be based on the following very simple 
theorem. 

Theorem 2. Let C/j and Ui, for i = 1, . . . , K , be independent families of cen- 
tered Gaussian random variables, whose covariances satisfy the inequalities 
for generic configurations 

EiUiUj) = S {j > E(U t Uj) = 4-, (25) 

and the equalities along the diagonal 

E(U t Ui) ee S u = E(UiUi) ee 4, (26) 

then for the quenched averages we have the inequality in the opposite sense 

ElogJ2wiexp(Ui) < ElogJ^WiexpiUi), (27) 

i i 

where the w $ > are the same in the two expressions. 

The proof is extremely simple and amounts to a straigthforward calcula- 
tion. In fact, let us consider the interpolating expression 

E\ogJ2 w ^ e MViU i + VT^Ui), (28) 

i 

where < t < 1 . Clearly the two expressions under comparison correspond to 
the values t = and t = 1 respectively. By taking the derivative with respect 
to t, and then integrating by parts with respect to the Gaussian variables, we 
immediately see that the interpolating function is nonincreasing in t, and the 
theorem follows. On the other hand, considerations of this kind are present 
in the mathematical literature of some years ago. Two typical references are 
[ISj and PH. 

We give here some striking applications of the basic comparison Theorem. 
In |Sj we have given a very simple proof of a long waited result, about the 
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convergence of the free energy per site in the thermodynamic limit. Let us 
show the argument. Let us consider a system of size N and two smaller 
systems of sizes Ni and N 2 respectively, with iV = N± + N 2 , as before in the 
ferromagnetic case. Let us now compare 

ElogZ N (P,h,J) = Elog exp((3^-K(aj)exp((3hY,^), (29) 

ai...U]\f i 

with 



(T1...(TJV 



= E log Z Nl (J3, h,J) + E log Z N2 ((3, h, J) , (30) 

where are the (cxj, i = 1, . . . , Ni), and are the (<7j, i = iV\-|-l, . . . , N). 
Covariances for and are expressed as in (|22|). but now the overlaps 
are substituted with the partial overlaps of the first and second block, qi and 
q 2 respectively. It is very simple to apply the comparison theorem. All one 
has to do is to observe that the obvious 

Nq = N iqi + N 2 q 2 , (31) 

analogous to (JOJ), implies, as in (JUJ), 

2 N x 2 N 2 2 
9 < -Jf9i + Jf<h- (32) 

Therefore, the comparison gives the superaddivity property, to be compared 
with ©, 

E log h,J)>E log Z M (/?, h,J) + E log (/?, /i, J) . (33) 
From the superaddivity property the existence of the limit follows in the form 
lim N^E log Z N (P, h, J) = sup N^E log Z N (f3, h, J), (34) 

AT— >oo jy 

to be compared with (|12)1. 

The second application is in the form of the Aizenman-Sims-Starr gen- 
eralized variational principle. Here, we will need to introduce some auxil- 
iary system. The denumerable configuration space is given by the values of 
a = 1,2,.... We introduce also a probability measure w a for the a system, 
and suitably defined overlaps between two generic configurations p(a,a'), 
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with p(a, a) — 1. A family of centered Gaussian random variables k(a), now 
indexed by the configurations a, will be defined by the covariances 

E(k(a)k(a')) = p 2 (a,a'). (35) 

We will need also a family of centered Gaussian random variables r)i(a), 
indexed by the sites i of our original system and the configurations a of the 
auxiliary system, so that 

E{r)i{a)rn>{a')) = 6 iV p(a,a'). (36) 

Both the probability measure w a , and the overlaps p(a, a') could depend 
on some additional external quenched noise, that does not appear explicitely 
in our notation. 

In the following, we will denote by E averages with respect to all random 
variables involved. 

In order to start the comparison argument, we will consider firstly the 
case where the two a and a systems are not coupled, so to appear factorized 
in the form 

<j\. ..(Tat ol i 

= E log Z N ((3, h,J) + E log ™* exp(/3^«(a)). (37) 

a 

In the second case the k fields are suppressed and the coupling between 
the two systems will be taken in a very simple form, by allowing the r\ field 
to act as an external field on the a system. In this way the a's appear as 
factorized, and the sums can be explicitely performed. The chosen form for 
the second term in the comparison is 

Elog ^2 ^ w<x exp(/3 ^ r}i(a)<Ti) exp((3h ^ Gj) 

ai...aj\f a i i 

= N\og2 + E\ogJ2 w a (ciC 2 ...c N ), (38) 

a 

where we have defined 

q = cosh j3(h + r/i (a)), (39) 
as arising from the sums over cr's. 
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Now we apply the comparison Theorem. In the first case, the covariances 
involve the sums of squares of overlaps 

~(q 2 (a,a')+p 2 (a,a')). (40) 

In the second case, a very simple calculation shows that the covariances 
involve the overlap products 

q(a,a')p(a,a'). (41) 

Therefore, the comparison is very easy and, by collecting all expressions, we 
end up with the useful estimate, as in ^2], holding for any auxiliary system 
as defined before, 

N~ X E log Z N ({3, h, J) < (42) 
log 2 + N~ 1 (E log ^ w a(cic 2 . . . c N ) - E log ^ w <* ex P(P]f^K a )))- 



5 The Parisi representation for the free en- 
ergy 

We refer to the original paper and to the extensive review given in jH], 
for the general motivations, and the derivation of the broken replica Ansatz, 
in the frame of the ingenious replica trick. Here we limit ourselves to a 
synthetic description of its general structure, independently from the replica 
trick 

First of all, let us introduce the convex space X of the functional order 
parameters x, as nondecreasing functions of the auxiliary variable q, both x 
and q taking values on the interval [0, 1], i.e. 

X3x: [0,1] 3q^x{q) € [0,1]. (43) 

Notice that we call x the function, and x(q) its values. We introduce a metric 
on X through the L x ([0, 1], dq) norm, where dq is the Lebesgue measure. 

For our purposes, we will consider the case of piecewise constant func- 
tional order parameters, characterized by an integer K, and two sequences 
go, qi, ■ ■ ■ , qx, mi, m.2, • • • , rriK of numbers satisfying 

= go < <?i < ■ ■ ■ < gx-i < qn — 1, < mi < m 2 < . . . < m K < 1, (44) 
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such that 



x(q) = mi for = q < q < q\, x(q) = m 2 for qi < q < q 2 , 

. . . , x(q) = m K for q K „ x < q < q K . (45) 

In the following, we will find convenient to define also m = 0, and m^+i = 1. 
The replica symmetric case of Sherrington and Kirkpatrick corresponds to 

K = 2, q\ = q, mi = 0, m 2 = 1. (46) 

Let us now introduce the function /, with values f(q, y\ x, j3), of the vari- 
ables q € [0,1], y G R, depending also on the functional order parameter x, 
and on the inverse temperature (3, defined as the solution of the nonlinear 
antiparabolic equation 

(dj)(q,y) + \(d 2 y f)(q,y) + l -x{q){d y f)\q,y) = 0, (47) 
with final condition 

/(l,y)=]ogcosh(/fy). (48) 

Here, we have stressed only the dependence of / on q and y. 

It is very simple to integrate Eq. (j^Tj) when x is piecewise constant. In 
fact, consider x(q) = m a , for q a _i < q < q a , firstly with m a > 0. Then, it is 
immediately seen that the correct solution of Eq. (jlTj) in this interval, with 
the right final boundary condition at q = q a , is given by 

/(<?, V) = — log / exp(m a f(q a , y + z^J q a - q)) dfi(z), (49) 
m a J 

where dfi(z) is the centered unit Gaussian measure on the real line. On the 
other hand, if m a = 0, then ()47|) loses the nonlinear part and the solution is 
given by 

f(q,y) = J f{q a ,y + z^q~^q)d^{z), (50) 

which can be seen also as deriving from (|4*9*|) in the limit m a — > 0. Starting 
from the last interval K, and using (j4*9*|) iteratively on each interval, we easily 
get the solution of (jUJ), (|4^j) . in the case of piecewise order parameter x, as 
in (1431) . 

Now we introduce the following important definitions. The trial auxiliary 
function, associated to a given mean field spin glass system, as described in 
Section 3, depending on the functional order parameter x, is defined as 

\og2 + f{0,h;x,P)-^^ qx{q)dq. (51) 
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Notice that in this expression the function / appears evaluated at q — 0, 
and y = h, where h is the value of the external magnetic field. This trial 
expression shoul be considered as the analog of that appearing in (fT3|) for 
the ferromagnetic case. 

The Parisi spontaneously broken replica symmetry expression for the free 
energy is given by the definition 

- Pfp{/3,h) =M(log2 + f(0,h',x,/3)-^ [ qx(q)dq), (52) 
x ^ Jo 

where the infimum is taken with respect to all functional order parameters 
x. Notice that the infimum appears here, as compared to the supremum in 
the ferromagnetic case. 

In jO], by exploiting a kind of generalized comparison argument, involving 
a suitably defined interpolation function, we have established the following 
important result. 

Theorem 3. For all values of the inverse temperature (3, and the external 
magnetic field h, and for any functional order parameter x, the following 
bound holds 

N~ l E log Z N (P, h, J) < log 2 + f(0, h; x, (3) - ^ j q x(q) dq, 
uniformly in N. Consequently, we have also 

N~ X E log Z N {(3, h, J) < inf (log 2 + f(0, h; x, 0) - ^ / q x(q) dq), 

x * Jo 

uniformly in N . 

However, this result can be understood also in the frame of the general- 
ized variational principle established by Aizenman-Sims-Starr and described 
before. 

In fact, one can easily show that there exist an a systems such that 
iV -1 ^ log ^ w a{cic 2 ■■■ c N ) = f(0, h; x, (3), 

a 

/W f3 2 f 1 

— k(a)) = — J qx(q) dq, 

uniformly in N. This result stems from previous work of Derrida, Ruelle, 
Neveu, Bolthausen, Sznitman, Aizenman, Talagrand, Bovier, and others, 
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and in a sense is implicit in the treatment given in We plan to deal with 
this important representation in a forthcoming note. 

We see that the estimate in Theorem El are also a consequence of the 
generalized variational principle. 

Up to this point we have seen how to obtain upper bounds. The problem 
arises whether, as in the ferromagnetic case, we can also get lower bounds, 
so to shrink the thermodynamic limit to the value given by the inf^ in The- 
orem 01 After a short announcement in ^Qj, Michel Talagrand wrote an 
extended paper jllj . to appear on Annals of Mathematics, where the com- 
plete proof of the control of the lower bound is firmly established. We refer 
to the original paper for the complete details of this remarkable achievement. 
About the methods, here we only recall that in jH] we have given also the 
corrections to the bounds appearing in Theorem El albeit in a quite compli- 
cated form. Talagrand, with great courage, has been able to establish that 
these corrections do in fact vanish in the thermodynamic limit. 

In conclusion, we can establish the following extension of Theorem ^ to 
spin glasses. 

Theorem 4. For the mean field spin glass model we have 

lim N~ X E log Z N (J3, h, J) = sup N" X E log Z N (J3, h, J) (53) 

N~>oo 

= inf(log2 + /(0,M,/3)-^ / qx{q)dq). (54) 
x * Jo 

6 Conclusion and outlook for future develop- 
ments 

As we have seen, in these last few years there has been an impressive progress 
in the understanding of the mathematical structure of spin glass models, 
mainly due to the systematic exploration of comparison and interpolation 
methods. However many important problems are still open. The most im- 
portant one is to establish rigorously the full hierarchical ultrametric or- 
ganization of the overlap distributions, as appears in Parisi theory, and to 
fully understand the decomposition in pure states of the glassy phase, at low 
temperatures. 

Moreover, is would be important to extend these methods to other impor- 
tant disordered models as for example neural networks. Here the difficulty 
is that the positivity arguments, so essential in comparison methods, do not 
seem to emerge naturally inside the structure of the theory. 

We plan to report on these problems in future works. 



15 



Acknowledgments 

We gratefully acknowledge useful conversations with Michael Aizenman, 
Pierluigi Contucci, Giorgio Parisi and Michel Talagrand. The strategy ex- 
plained in this paper grew out from a systematic exploration of compari- 
son and interpolation methods, developed in collaboration with Fabio Lucio 
Toninelli. 

This work was supported in part by MIUR (Italian Minister of Instruc- 
tion, University and Research), and by INFN (Italian National Institute for 
Nuclear Physics). 

References 

[I] D. Sherrington and S. Kirkpatrick, Solvable Model of a Spin-Glass Phys. 
Rev. Lett. 35, 1792-1796 (1975). 

[2] S. Kirkpatrick and D. Sherrington, Infinite-ranged models of spin-glasses, 
Phys. Rev. B17, 4384-4403 (1978). 

[3] D.L. Stein, Disordered Systems: Mostly Spin Glasses, in: Lectures in the 
Sciences of Complexity, ed. D.L. Stein, Addison- Wesley, NY, 1989. 

[4] M. Mezard, G. Parisi and R. Zecchina, Analytic and Algorithmic Solution 
of Random Satisfiability Problems, Science 297, 812 (2002). 

[5] M. Mezard, G. Parisi and M. A. Virasoro, Spin glass theory and beyond, 
World Scientific, Singapore, 1987. 

[6] M. Talagrand, Spin glasses: a challenge for mathematicians. Mean field 
models and cavity method, Springer- Verlag, Berlin (2003). 

[7] F. Guerra, Sum rules for the free energy in the mean field spin glass 
model, Fields Institute Communications 30, 161 (2001). 

[8] F. Guerra and F. L. Toninelli, The Thermodynamic Limit in Mean Field 
Spin Glass Models, Commun. Math. Phys. 230, 71-79 (2002). 

[9] F. Guerra, Broken Replica Symmetry Bounds in the Mean Field Spin 
Glass Model, Commun. Math. Phys. 233, 1-12 (2003). 

[10] M. Talagrand, The Generalized Parisi Formula, Compte Rendu de 
lAcademie des Sciences, Paris 337, 111-114 (2003). 

[II] M. Talagrand, The Parisi formula, Annals of Mathematics, to appear. 



16 



[12] M. Aizenman, R. Sims and S. Starr, Extended variational principle for 
the Sherrington- Kirkpatrick spin-glass model, Phys. Rev. B68, 214403 
(2003). 

[13] H.E. Stanley, Introduction to phase transitions and critical phenomena, 
Oxford University Press, New York and London, 1971. 

[14] D. Ruelle, Statistical mechanics. Rigorous results, W.A. Benjamin Inc., 
New York, 1969. 

[15] K. Joag-dev, M.D. Perlman and L.D. Pitt, Association of normal ran- 
dom variables and Slepian's inequality, Annals of Probability 11, 451-455 
(1983). 

[16] J.-P. Kahane, Une inegalite du type Slepian and Gordon sur les processus 
gaussiens, Israel J. Math. 55, 109-110 (1986). 

[17] G. Parisi, A sequence of approximate solutions to the S-K model for 
spin glasses, J. Phys. A13, L-115 (1980). 



17 



